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ABSTRACT 

1. earning  curves  have  been  used  extensively  for  predictive  pur¬ 
poses  in  the  airframe  and  other  industries.  In  many  instances  this  has 
led  to  erroneous  results  because  analysts  failed  to  extend  learning 
curve  theory  and  dev-lop  adequate  analytical  techniques  in  the  turbulent 
regime  of  the  cost  history  characterizing  these  industries.  It  is  this 
area  where  a  series  of  design  changes  induces  a  seriesof  perturbations 
whose  turbulence  intensity  is  a  function  of  iiic-  fr  quency  of  occurrence 
and  magnitudes  of  the  design  changes  under  consideration. 

so  l\V<  — •_  { / 

-In  1  T,  a  series  of  formulations  amenable  to  machine  program¬ 
ming  was  developed  for  the  accurate  determination  of  perturbed  unit 
costs.  This  development  was  based  on  additions  of  new  learning  having 
a  constant  slope.  .  .  .  „ . 

i 

In  this  discussion,  the  development  of  ftefr  i-will  be  generalized  by 
developing  formulas  for  the  addition  of  new  effort  having  variable 
slopes.  Consideration  will  also  be  given  to  the  expressions  involving 
elementary  unit  cost  expressions  so  that  cumulative  average  and  cumu¬ 
lative  total  values  can  readily  be  obtained  from  existing  experience 
curve  tables.  Conversely,  the  problem  of  determining  the  magnitudes 
of  design  changes  and  the  slopes  of  new  effort  from  graphical  data  will 
also  be  considered.  ( 

■V 

INTRODUCTION 

The  purpose  of  this  investigation  is  to  generalize  and  extend  the  expressions  developed 
in  Ref.  1.  In  the  interest  of  generalization,  formulas  for  the  addition  of  new  effort  having  vari¬ 
able  slopes  will  be  developed.  Conversely,  a  method  will  be  presented  for  the  confirmation  of 
the  slope  of  the  additional  new  learning  and  the  magnitudes  of  the  design  changes  imposed.  In 
addition,  the  analytical  expression  for  the*  1th  perturbation  will  be  expressed  in  terms  of  the  i 
design  change  parameters  and  i-t- 1  elementary  functions.  The  latter  consist  of  the  initial  unit 
cost  expression  and  the  i  unit  cost  expressions  resulting  from  the  additions  of  new  effort  as¬ 
sociated  with  the  l  design  changes  which  follow.  With  this  expression  for  the  i**1  perturbation 
in  terms  of  (i  +  1)  unit  cost  expressions,  it  will  be  possible  to  develop  cumulative  average  and 
cumulative  total  cost  evaluations  based  on  experience  curve  tables  for  use  at  any  point  in  the 
turbulent  regime. 

DERIVATION  OF  APPLICABLE  EQUATIONS 

The  choice  of  symbolism  in  Ref.  1  is  inconsistent  with  that  in  common  usage.  The  first 
requirement,  therefore,  Is  to  re-define  quantities  so  that  this  inconsistency  will  be  less  pro¬ 
nounced.  As  a  start,  the  Initial  unit  cost  expression  will  be  changed  to 

595 


BEST  AVAILABLE  COPY 


.39  t'. 

(1) 

A. “.erf 

Ug  =  cost  of  first  production  item 
bg  =  numerical  value  of  exponent 
x  =  production  quantity. 

The  subscript  zero  in  the  foregoing  characterizes  cost  conditions  in  the  unperturbed  regime. 

The  turbulent  cost  regime,  on  the  other  hand,  is  characterized  by  a  series  of  cost  dis¬ 
continuities  induced  by  a  series  of  major  design  changes.  The  degree  of  turbulence  is  a  func¬ 
tion  of  the  frequency  of  occurrence  of  these  design  changes  and  their  magnitudes. 

The  formulas  associated  with  the  turbulent  cost  regime  are  characterized  by  the  sub¬ 
script  1  where  i=  1,2,  .  .  . ,  n,  when  evaluated,  is  indicative  of  a  particular  perturbation  caused 
by  a  particular  design  change.  The  ith  perturbation  will  become  effective  after  the  completion 
of  item  X(  and  continue  in  effect  until  item  X.  +  ^  is  completed.  Associated  with  the  i'*1  major 
design  change  are  two  weighting  factors  A^  and  t.  along  with  the  new  parameter,  b;  ,  where 
the  latter  is  indicative  of  variable  new  learning  capabilities,  a^  represents  the  present  reduc¬ 
tion  in  the  previous  manufacturing  effort,  f.  jCX-Xj  j).  while  t  represents  the  percent  of 
new  effort  to  be  added.  The  unit  cost  expression,  T(X-X.),  after  the  ith  discontinuity  in  the 
range,  Xj  <  X  --  Xi  +  1.  is  given  by  the  following  recursion  formula: 

(2)  f.(X-Xj)  =  (1  - Aj) -Xj.j)  +  aQ^X -XjJ  1 

where  (i  =  1,  2 . n)  and  Xn  =  0. 

v  l  L 

it  is  now  possible  to  state  the  expression  for  the  il  perturbation  in  terms  of  the  i  de¬ 
sign  change  parameters  and  i+  1  elementary  functions.  As  indicated  previously  the  latter  con¬ 
sist  of  the  initial  unit  cost  expression  and  the  i  unit  cost  expressions  resulting  from  the  addi¬ 
tions  of  variable  new  effort  associated  with  the  i  design  changes  which  follow.  By  using  Eqs. 
(1)  and  (2)  and  substituting  f^jfX-X^j)  into  fj(X-Xj).  it  is  apparent  that  for  i  =  1.  2.  3.  etc., 
we  have: 

/  \“bj 

(3)  fj(X-Xj)  -  ( 1  -  A  j )  fQ(X)  *  Tj  agfx-X,] 

-  U-AjjagX'V  Tj  a0(X-Xlfl 


P.  \1  .IAMI-.S 

>0<x)  =  a0X  b° 


(4)  f2(X-X2)  -  (l-A2)f1(X-X1)  *  '  2  a0  X  -  x2 


=  (1-a2)  AX-^)^x'b°  .  h  a0  X-X1'blrT2a0i  X-X2rb2 


-b_  ,-b.  ,  ,-b_ 

*  (»- VT!  a0  X'Xl!  *  T2  a0(X  *  X2t 


*  (J.*2)(l-*,)a0X 
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(5) 


f3(X-X3)  =  (1  -  A3)f2(X-X2)  +  r3  a0(X-X3, 


*b3 


”fc>n 


(1  -  A3)  |(1  -  A2)(l-  AjJagX  u  +  (1  -  a2)  Tj  aQ ; X  -  X j 


-bo 


+  T2a0\X‘X2j  +  7 3  a0lA'A3.i 


,  Y  Y  \  3 

in  X-XJ 


-bn  /  'rbi 

=  (1-  A3)(l-  A2)(l  -  Aj)a0X  b  +  (1- A3)(l- A2)7j  a0'X-Xj  1 


~°2  /  '-t>3 

(1-A3)  r2a0(X-X2)  /  +  73a0(x-X3J 


Equation  (5)  can  be  expressed  more  compactly  as 

(6)  f3(X-X3)  =  C3(x-XQ'  b°  +  C2(x-xJ  bl  *  CuX-X2)b2  ♦  t3  aQ(x-X3)  ^ 

Equation  (6)  is  a  special  case  of  the  following  expression: 

J-M 


(7) 

where 

(8) 

and 

(9) 


-b. 


MX-X1>  =  2_,  ci-)  (x-xjP*  Tiao\x-xi 


rb; 


j=0 


Ci-i  =  a»  T‘ 


k=i 

o  ^  n v- 

k=j+l 


=  1. 


By  substitution  of  Eq.  (8)  in  Eq.  (7),  it  is  seen  that 

)=!-!  k=i 


(10) 


and 


(11) 


•.«  -  v  ■  *o  21  ’>  n  °  •  -  xifbi  •  *o  *  -  v** 


)=0  k= j  + 1 

j=i-l 


k=i 


fi(x-xi>  -  ao  y 7,  x-x,’bi  7T(,-*k,t  r. 

k=j  +  l 


a„;X-X 


\“i 


i  ‘or  i; 


In  connection  with  the  foregoing  development,  attention  is  invited  to  the  fact  that  either 
Eq.  (7)  in  conjunction  with  Eq.  (8)  or  their  composite  equivalent,  Eq.  (11).  is  an  expression  for 


598 


P  M  TAMPS 


the  ith  perturbation  in  terms  of  the  initial  unit  cost  expression,  Eq.  (1),  and  the  i  unit  cost 
expressions  associated  with  the  i  design  changes  which  follow.  The  latter,  it  will  be  noted, 
become  effective  witl.  production  quantities  Xj  ♦  !  where  i=l,2,  ...,n  respectively.  Equations 
(5).  (4)  and  (3)  are  evaluations  of  Eq.  (11)  or  Eqs.  (7)  and  (8)  for  i  =  3,2,  1  respectively.  A  veri¬ 
fication  of  Eq.  (11)  is  given  in  Appendix  A. 

Since  Eqs.  (7)  and  (8)  or  the  composite  relationship,  Eq.  (11),  involve  (i  +  1)  unit  cost 
expressions,  the  corresponding  cumulative  average  and  cumulative  total  values  can  readily  be 
obtained  from  existing  experience  curve  tables.  The  latter  are  particularly  useful  in  arriving 
at  an  equitable  financial  adjustment  for  the  extension  of  a  contract  and  the  procurement  of  addi¬ 
tional  production  quantities  in  the  turbulent  regime. 

DETERMINATION  OF  NEW  EFFORT  SLOPE  AND  DESIGN  CHANGE 
PARAMETERS  FROM  GRAPHICAL  DATA 

The  problem  of  determining  the  magnitudes  of  the  slope  and  design  change  parameters 
from  graphical  data  will  now  be  discussed.  A  requirement  for  the  determination  is  that  these 
quantities  must  be  known  for  each  of  the  preceding  i  perturbations  along  with  the  slope  of  the 
initial  unit  cost  expression.  In  cases  where  this  information  is  not  available,  it  will  be  neces¬ 
sary  to  start  with  the  graph  of  the  initial  unit  cost  expression  and  determine  the  design  change 
parameters  and  new  learning  slope  associated  with  the  first  perturbation.  With  this  informa¬ 
tion,  the  same  procedure  can  be  applied  to  the  2nd,  3rd,  and  finally  to  the  (i  +  1st)  perturbation. 
The  previously  derived  expressions  can  then  be  used  in  conjunction  with  a  technique  of  succes¬ 
sive  approximations. 

As  an  illustration  of  the  foregoing,  a  method  of  determining  the  design  change  parame¬ 
ters  will  be  illustrated  by  considering  for  simplicity  a  single  perturbation.  It  is  assumed  that 
the  initial  cost  of  an  item,  aQ,  is  100  and  that  bQ  =  0.322  for  801  learning.  Alter  the  comple¬ 
tion  of  production  item  X  =  30,  a  design  change  calling  for  ij  =  0.113  and  Tj  =  0.207  becomes 
effective.  The  value  of  bj  is  assumed  to  be  0.454  for  73%  learning. 

The  foregoing  design  change  is  represented  graphically  on  the  linear  scale,  Figure  1, 
and  the  conventional  log-log  scale,  Figure  2.  Without  being  influenced  by  tnese  design  changes, 
the  problem  now  is  to  use  graphical  values,  and  compute  values  of  >j  .  Tj  and  b.  .  With  this 
in  mind,  the  coordinates  of  three  points  read  from  the  graph  had  the  following  values: 


X 

fj(X-X, 

31 

50.0 

32 

44.2 

34 

39.8 

When  these  three  -,e»s  of  values  are  substituted  in  the  governing  equation 

-*>n  '  rhj 

1  j(X  - X)  -  (1  -  *,)aQX  u  *  aQ  Tj  X-Xj 

three  equations  in  three  unknowns  are  obtained.  Because  of  the  transcendental  nature  of  this 
set  of  equations  together  with  the  round-off  errors  and  errors  resulting  from  using  graphical 
values,  the  results  ol  the  computations  given  in  Appendix  B  left  something  to  be  desired. 
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In  an  effort  to  obtain  acceptable  accuracies  in  the  computation  of  Aj  ,  Tj  ,  and  bj  ,  it 
was  necessary  to  resort  to  a  method  of  successive  approximations.  In  this  connection,  the  ap¬ 
plicable  expression  fj(X  -Xj)  is  unitized  and  rearranged  so  that  the  new  effort  is  expressed  as 
a  function  of  the  graphical  value,  the  initial  unit  cost  expression,  and  the  parameter  >j  .  Points 
are  read  from  the  graph  in  the  neighborhood  of  the  discontinuity  and  one,  a  check  point,  is 

taken  as  far  from  the  discontinuity  as  possible.  B>  inserting  trial  values  of  A.  =  0.05,  0.10, 

-b, 

015,  etc.,  re.  responding  slopes,  S  =  2  1 ,  are  obtained.  On  the  basis  of  the  slope  which  pro¬ 

duces  a  prediction  whose  deviation  from  the  extreme  check  point  is  a  minimum,  values  of  Aj 
and  Tj  are  computed.  By  using  values  of  Aj,  in  the  neighborhood  of  the  value  just  computed, 
the  program  is  reptated  until  convergence  develops.  Details  of  the  method  follow: 

-bn  -h. 

(12)  fj(X-Xj)  =  (1-  AjUqX  +  aQ  Tj  ;X  -  X j  \ 

(13)  gjtX-Xj)  =  f^X-Xj)  aQ  =  (1-  A,)X  b°  *  ^(x-Xjj  ^  , 
and 

(14)  Tl(X‘xl'  1  =  g^X-Xj)  -  (1-  Aj)x  °, 

-t>  i 

where  7j  X  -Xj  is  the  new  effort  in  terms  of  the  graphical  value,  gj(X  -Xj),  the  original 

unitized  unit  cost  expression,  X  and  the  parameter  A^ .  By  using  values  X,  X",  X"'  where 
X’  =  Xj  +  1.  X"  =  Xj  +  2,  X"  1  =  Xj  +  4,  Eq.  (14)  produces  *he  following: 

(15)  g^X’-Xj)  -  (1-  Ajjfx1:  b°  =  Tj  , 

(16)  gj(X"-Xj)  -  (1  -  Aj)(x”)  b°  =■  Tj (2)  b*  , 

and  -b  -b 

(17)  g1(X"'-X1)  -  (1-  Aj)(x"’)  °=  7,(4)  1  . 

Equations  (15),  (16),  and  (17)  are  expressed  in  tabular  form  for  Aj  =  0.05,  0.10,  0.15,  respec¬ 
tively,  in  the  following  Table  1: 


TABLE  1 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(2)  -  (4) 

(8) 

(2)  -  (5) 

(9) 

(2)  -  (6) 

X 

gjtX-X,) 

x'J° 

'b0 
0.95X  u 

-bo 
0.90X  u 

”b0 
0.85X  u 

/  rbi 

Ai  ix  ■  xi) 

rbi 

AiVx-xJ 

Al(X-Xl)’bl 

Aj  =  0.05 

Aj  =  0.10 

>,  =  0.15 

31 

0.500 

0.3311 

0.2980 

0.2814 

WM 

0.2020 

0.2186 

32 

0.442 

0.3277 

0.3177 

0.2949 

0.2785 

I 

0.1471 

34 

0.398 

0.3053 

0.2893 

0  2732 

■SB 

0.1087 

0.345 

0.2858 

0.2745 

0.2592 

0.0705 

D 

0.265 

0.2547 

0.2420 

0.2292 

0.2165 

Ea 

0.0358 

■ 
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Division  of  Eq.  (16)  by  (15)  produces  the  slope  S  =  2  The  slope  numerical  values  for 
Aj  =  0.05,  0.10,  0.15,  respecti-ely,  are  obtained  by  dividing  the  second  line  by  the  first  in  col¬ 
umns  (7),  (8)  and  (9)  of  Table  1,  therefore,  Table  2. 


TABLE  2 
(PARTS  1  AND  2) 


PART  1 

(10) 

(ID 

(12) 

(13) 

(14) 

1 

f2  *bi 

S  =  —  =  2 
<1 

0.1855(x-Xj)  bl 
Aj  =  0.05 

0.2020  (X-Xj  1 

A.  =  0.10 

0.2186[X-Xj)  1 
Aj  =  0.15 

Comments 

Xt  =  0.30;  X  =  40 

Xj  =  0.30;  X  =  40 

Xj  =  0.30;  X  =  40 

0.05 

0.7046 

0.0567 

Columns  (12) 

(13)  and  (14)are 

0.10 

0.7272 

0.0710 

predictions 

based  on 

0.15 

0.7479 

_ 

__ 

0.0841 

slopes,  col.  (11) 
for  Aj  =  0.05, 
0.10,  and  0.15 

PART  2 

0.05 

0.0552 

These  values 

tvere  taken  from 

0.10 

0.0705 

Cols.  (7),  (8), 
and  (9)  with 

0.15 

0.0858 

o 

II 

X 

_ 

0.0015 

0.0005 

-0.0017 

In  the  first  section  of  the  foregoing  table  covering  columnB  (10)  io  (14)  predictions  are 
made  for  the  cost  of  the  new  effort  at  X  =  40  for  slopes  of  0.7046,  0.7282  and  0.7479  based  on 
initial  trial  values  of  Aj  =  0.05,  0.10,  0.15  respectively.  The  coefficients  of  the  negative  ex¬ 
ponential  In  the  new  effort  expressions,  columns  (12),  (13)  and  (14)  are  the  numbers  given  on 
the  first  line  under  columns  (7),  (8)  and  (9).  The  values  appearing  in  the  second  section  of  this 
table  under  columns  (12),  (13)  and  (14)  are  the  values  for  X  -  40.  line  4.  under  columns  (7), 

(8)  and  (9).  From  the  differences  given  in  the  last  line,  it  is  apparent  that  the  slope  lies  be¬ 
tween  0.7282  and  0.7497.  The  value  0.74  will  be  used  for  the  second  Iteration. 

By  substituting  values  for  X  =  31  and  32  under  columns  (2)  and  (3)  into  Eqs.  (15)  aru 
(16),  we  obtain 
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(18) 

0.500 

-  0.33 1 1  (x  -  Aj)  =  r 

and 

(19) 

0.442 

-  0.3277(1  -  A  )  =  (1.74)  t. 

with  the  solution  being  given  by  =  0.2118  and  A^  =  0.1294. 

The  second  iteration  was  based  on  Aj  values  of  0.10,  0.13  and  0.16.  These  were  used 
in  columns  (4),  (5)  an  1  (6)  respectively  and  the  entire  computation  repeated.  At  the  end  of  the 

-b, 

second  iteration,  the  following  numerical  values  were  obtained:  S  =  2  -  0.736,  =  0.2083 

and  Aj  =  0.121. 

Ttv  third  iteration  was  based  on  a^  values  in  the  neighborhood  of  0.121  which  was  the 

result  ol  the  second  iteration.  The  actual  a  values  substituted  in  columns  (4),  (5)  and  (6)  were 

-b, 

0.11,  0.12  and  0.13.  This  iteration  produced  the  following  values:  S  =  2  1  =  0.732,  t  =  0.2053 

and  *.  (.average)  =  C.112.  The  latter,  it  will  be  observed,  compares  very  lavorablv  with  the 
*bl 

true  values  of  S  =  2  1  0.73,  =  0.207  and  Aj  =  0.113. 

RESULTS 

The  iormulations  derived  herein  are  suitable  for  accurately  describing  the  entire  his¬ 
tory  of  l  production  item  in  the  turbulent  regime  when  additions  of  new  learning  are  individu¬ 
ally  different:  that  is.  the  formulations  account  for  i  different  learning  rates  when  i  design 
changes  are  specified.  These  equations  have  also  been  expressed  in  terms  of  (i  *  1)  elemen¬ 
tary  unit  cost  expressions  so  tiiat  cumulative  average  and  cumulative  total  values  can  readily 
be  obtained  from  existing  tables  for  contract  adjustments.  By  using  a  technique  of  successive 
approximati  ms.  the  exp  est  ions  can  be  used  in  conjunction  with  graphical  data  for  the  deter¬ 
mination  of  not  only  the  design  change  parameters  but  the  learning  rate  associated  with  the 
addition  of  new  effort  when  a  design  change  has  been  made.  All  formulations  are  well  suited  to 
machine  programming  ciforts  and  repetitive  computational  procedures. 


Ap^’ndix  A 

VERIFICATION  OF  EQUATION  (11) 


l-t-i 


k-i 


f,(X  -  X 


•  •  X:  :":"77 


d-V+  Vo\x*xt 


fbi 


1=0 


k= j  *  1 


Letting  i  3 .  it  is  seen  that 


-b  ^-3  -b 

f3(X-x3>  ,  a0  <;VX-X0  0  77  Tl(x  *XI'  1  77<l-V 

!  k= 1  k=2 

t-b  1  -t) 

*  -2  X-X2)  2  77  0- v.v*  ’3  a0  X  ’  X3  3 

k=3 


ana:  yv.s  or  proc.kkss  .  ri<'-K 


GC3 


(1A;  f3(X-X3)  =  ao  TOlX-Xo'’b°(|-Al)(,-A2)(1->31 

,-b, 

♦  a0  *i  \X  *  xi;  -  V  (1-  a3> 

/  \'bo 

+  aQ  7 2  lx  "  x2;  (]  ‘  a3> 

i  r  b  q 

*  a0  t3  tX  "  x3 

In  View  o{  the  fact  that  :Q  =  1  by  Eq.  (9).  it  is  apparent  that  Eq.  (1A)  is  the  same  a 

Eq.  (U). 


Appendix  B 

COMPUTATIONS  USING  GRAPHICAL  VALUES 


(IB) 

f,(X-X.) 

=  (1-  A,)a0X  VaQ  yX-X 

(2B) 

(t.'X-X.) 

A.  1 

=  f,(X-Xj)  ,  (1-  A()x‘b° 

X 

fjOC-X.) 

fl<X-X,)/a0 

O 

1 

X 

31 

50.0 

0.500 

0.3310 

32 

44.2 

0.442 

0.3277 

34 

39.8 

0.398 

0.3214 

Substitution  o!  the  above  tabular  values  in  Eq.  (2B)  produces. 

0.5000  =  0.3310(1  -  Aj)  *  7lj'l')*bl 

0.4420  =  0.3277(1  -  A,)  +  t,  2  ~ ^ 1 

,-b . 

0.3S80  =  0.3214(1  -  Aj)  ♦  r  4  1 

and 

0.1690  =  -0.3310  a  .  7j 

-b. 

0.1143  -■  -0.3277  a  *  Tj  2  1 

•b . 

0.0766  =  -0.3214  A  *  •  4  5 


